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Abstract. We report our Monte Carlo results on the critical and multicritical behavior of the 
±J Ising model [with a random-exchange probability P(J xy ) = pS(J xy — J) + (1 — p)S(J xy + J)], 
in two and three dimensions. We study the transition line between the paramagnetic and 
ferromagnetic phase, which extends from p = 1 to a multicritical (Nishimori) point. By a finite- 
size scaling analysis, we provide strong numerical evidence that in three dimensions the critical 
behavior along this line belongs to the same universality class as that of the critical transition 
in the randomly dilute Ising model. In two dimensions we confirm that the critical behavior is 
controlled by the pure Ising fixed point and that disorder is marginally irrelevant, giving rise 
to universal logarithmic corrections. In both two and three dimensions, we also determine the 
location of the multicritical Nishimori point, as well as the renormalization-group dimensions 
of the operators that control the renormalization-group flow close to it. 



1. Introduction 

The ± J Ising model is a simplified model [1] for disordered uniaxial magnetic systems showing 
glassy behavior in some region of their phase diagram and represents an interesting theoretical 
laboratory to study the effects of quenched disorder and frustration. It is defined by the lattice 
Hamiltonian 

H = — JxyCTxVy, (1) 

(xy) 

where a x = ±1, the sum is over the nearest-neighbor sites of a simple cubic lattice, and the 
exchange interactions J xy are uncorrelated quenched random variables, taking values ±J with 
probability distribution 

P(Jxy) = pS(J xy - J) + (1 - P )S(J xy + J). (2) 

In the following we set J = 1 without loss of generality. For p = 1, we recover the standard 
Ising model, while for p = 1/2, we obtain the usual bimodal Ising spin-glass model. The phase 
diagram in two and three dimensions is sketched in Figure [TJ it is symmetric for p — > 1 — p and 
thus here and in the following we only consider the case 1 — p < 1/2, i.e., p > 1/2. 




Figure 1. Phase diagram of the ±J Ising model in the T-p plane, in dimensions d = 2 (left) 
and d = 3 (right). 

In both d = 2 and d = 3, the model exhibits a paramagnetic- ferromagnetic (PF) phase 
transition in the region of low frustration. The PF transition line starts at the Ising point 
Xi s = (T = Ti s ,p = 1), where T\ s is the critical temperature of the Ising model, and extends 
up to the multicritical Nishimori point (MNP) at Xmnp = (T*,p*). Along this line, the critical 
behavior is analogous to that observed in randomly dilute Ising (RDI) models. 

In two dimensions, renormalization-group (RG) and conformal field theory predict j2j[3jlU[5] 
that disorder is marginally irrelevant, thus the critical behavior is controlled by the pure 
Ising fixed point, with universal logarithmic corrections [6]. In three dimensions disorder is 
relevant and the model shows a second-order phase transition in the three-dimensional (3D) 
RDI universality class [7], which describes transitions in generic diluted Ising systems with 
ferromagnetic exchange interactions (Sj. 

As argued in [U [THl E] , the MNP is located along the so-called Nishimori line (N line) [321 [13] 
defined by the equation 

tanh/? = 2p-l, (3) 

where (5 = 1/T. At the MNP the transition line is predicted to be parallel to the T axis |llj . 
Then, it reaches the T = axis at X c = (0,p c ). It has been proved [12] that ferromagnetism can 
only exist in the region p > p*; thus, p c must satisfy the inequality p c > p*. Recent numerical 
works indicate that p c is strictly larger than p*, though deviations are quite small, both in two 
|14j and three [15] dimensions, as well as in related models \16\ I17j. 

The critical behavior for p < p* depends on the dimension. The 3D ±J model exhibits a 
paramagnetic-glassy (PG) transition line, which extends from the MNP up to p = 1/2. On 
this transition line, the critical behavior is independent on p and belongs to the Ising spin-glass 
universality class [18J. In contrast, in two dimensions there is no evidence of a finite-temperature 
glassy phase. Glassy behavior is only expected for T = and p < p c : the glassy phase at T = 
is unstable with respect to thermal fluctuations. 

2. Finite-size scaling 

2.1. Paramagnetic- ferromagnetic transition line 

We present here the strategy used to analyze the critical behavior on the PF transition line in 
the 3D ± J model. The two-dimensional (2D) case is analyzed in [6]. 

According to the renormalization group (RG), in the case of periodic boundary conditions 
and for L — ► oo, where L is the lattice size, a generic RG invariant quantity R at the critical 



temperature 1//3 C behaves as 

R(L, (3 = /3 c ) = R* (1 + c u L- w + c 12 L~ 2uj + ■■■ + ^L'" 2 + •••), (4) 

where R* is the universal infinite- volume limit and u and oo 2 ar e the leading and next-to-leading 
correction-to-scaling exponents. In 3D RDI systems scaling corrections play an important role, 
since uj is quite small; indeed, we have u = 0.29(2) and oj 2 = 0.82(8) |191 l8]. Instead of computing 
the various quantities at fixed Hamiltonian parameters, we keep a RG invariant quantity R fixed 
at a given value Rf [20]. This means that, for each L, we determine the pseudocritical inverse 
temperature f3f(L) such that R{j3 = f3f(L),L) = Rf. All interesting thermodynamic quantities 
are then computed at = (3f(L). The pseudocritical inverse temperature (3f(L) converges to (3 C 
as L — > 00. The value Rf can be specified at will, as long as Rf is taken between the high- and 
low-temperature fixed-point values of R. The choice Rf = R* (where R* is the critical-point 
value) improves the convergence of f3f to [3 C for L — > 00; indeed [3f — [5 C = 0(L~ l l u ) for generic 
values of R f , while /3 f -/3 c = 0(L^ l / v ^) for R f = R* . 

We can then consider any other RG invariant quantity R a at fixed R = Rf, i.e., R a (L) = 
R a (L, fif(L)). For L — > 00, one can show that [8] 

R a (L) =R*(1 + b lx IT" + b 12 L~ 2 " + ■■■ + b 21 L'^ + •••), (5) 

R' a (L) = aoL 1 ^ (1 + a X iL- w + a 12 L- 2uj + ■■■ + a 21 L'^ + •••), (6) 

X(L) = X (L, [3 = (3 f (L)) = d L 2 -i (l + d lx L^ + d l2 L- 2w + ■■■ + d 2l L~^ + •■■)+ d b , (7) 

where R' a is the derivative of R a with respect to j3 and x is the susceptibility at fixed R = Rf. 
More details on this method can be found in [20\ [8]. 

2.2. Multicritical Nishimori point 

In the absence of external fields, the critical behavior at the MNP is characterized by two relevant 
RG operators. The singular part of the disorder-averaged free energy in a volume L d can be 
written as 

F sing (T,p,L) = L- d f( Ul Ly\u 2 L^,{ Ui L^}), i > 3, (8) 

where y\ > y 2 > 0, yi < for i > 3, Ui are the corresponding scaling fields, u\ = u 2 = at 
the MNP. The scaling fields U{ are analytic functions of the model parameters T and p. Using 
symmetry arguments, [10} [TT] showed that the scaling axis corresponding to u 2 = is along 
the N line, so that u 2 = tanh/? — 2p + 1. As for the scaling axis u\ = 0, e = 6 — d expansion 
calculations predict it [TT] to be parallel to the T axis. The extension of this result to lower 
dimensions suggests u\ = p — p* . Note that, if this conjecture holds, only the scaling field u 2 
depends on the temperature T. 

These results give rise to the following predictions for the FSS behavior around T*, p*. Along 
the N line, an RG invariant quantity R has the following behavior for L — > 00: 

R N = R* + b llUl L y i + ..., (9) 

where the subscript iV indicates that R is restricted to the ./V line and we have neglected scaling 
corrections. Its derivative R' with respect to (3 behaves as 

R' = bnu'^ 1 + b 2 W 2 L m + ■■■ = b 21 u' 2 Ly 2 + ..., (10) 

where the last equivalence holds only if u\ does not depend on the temperature. This result 
gives us a method to verify the conjecture of [11] : once y\ has been determined from the scaling 
behavior of a RG invariant quantity R close to the MNP, it is enough to check the scaling 
behavior of R' . If R' scales as L*> with £ < yi, the conjecture is confirmed and y 2 = £• Finally, 
along the ./V line the susceptibility x is expected to behave as 

XN = d L 2 -^ (1 + d x u x m + • • • ) . (11) 



3. Results 

3.1. The 3D ±J model at the paramagnetic-to-ferromagnetic transition line 

In [7] we perform MC simulations of Hamiltonian (pQ) at d = 3 close to the PF transition line. 
Using the method outlined in Sec. 12.11 we fit the large- L limit of various RG-invariant quantities 
at fixed R^ = ^/L = 0.5943, where £ is the correlation length; this value is very close to the fixed- 
point value i?| = 0.5944(7) [8] at (3 C - Our FSS analysis provides strong evidence that the critical 
behavior of the 3D ± J Ising model along the PF line belongs to the 3D RDI universality class. 
Indeed, all the RG-invariant quantities considered are in agreement with the results reported 
in [8]. Moreover, we find u = 0.682(3) and n = 0.036(2), in good agreement with the presently 
most accurate estimates [8] v = 0.683(2) and r\ = 0.036(1) for the 3D RDI universality class. 

We also note that the random-exchange interaction in the =LJ Ising model gives rise to 
frustration, while the RDI universality class describes transitions in generic diluted Ising systems 
with ferromagnetic exchange interactions. Therefore, our results imply that frustration is 
irrelevant along the PF transition line. Moreover, the observed scaling corrections are consistent 
with the RDI correction-to-scaling exponents uj = 0.29(2) and L02 = 0.82(8) [8]. Thus, 
frustration does not introduce new irrelevant perturbations with RG dimension \yf\ < 1. 

3.2. Multicritical Nishimori point 

In [141 115] . we perform MC simulations of Hamiltonian ([!]) at d = 2, 3, along the ./V line defined 
by ©. In both cases, the position of the MNP as well as the RG dimension y\ of the leading 
relevant operator u\ are determined by fitting several RG-invariant quantities R to ([9j) (in the 
2D case we also consider scaling corrections). We verify the conjecture that u\ does not depend 
on T and compute the RG dimension associated with u\ and U2- 

For the 2D model, the MNP is located at p* = 0.89081(7), T* = 0.9528(4), the RG dimensions 
of the operators that control the multicritical behavior are y\ = 0.655(15) and 2/2 = 0.250(2), 
and the susceptibility exponent is 77 = 0.180(5). In three dimensions, we find p* = 0.76820(4), 
T* = 1.6692(3), yi = 1.02(5), y 2 = 0.61(2), n = -0.114(3). 
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